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ABSTRACT 


A proof-of-concept experiment for remote sensing of precipitation by SAR is part of 
the SIR-C/X-SAR experiment. This thesis presents a feasibility study and 
recommendations for detection of precipitation using SIR-C/X-SAR. The principal 
limitation to rain measurement from a spacebome SAR is the poor SCR (signal-to- 
clutter ratio). This is in part due to the system configuration and largely due to the 
large magnitude of echoes associated with the surface component. Two geometries 
apply; off-vertical and vertical pointing angles. Here we present calculations for both. 

With vertical geometry a large clutter component is associated with range 
sidelobes of the chirped transmitter pulse. To overcome this problem a narrow 
transmitted pulse (3 psec) processed without dechirping was used. Since the 
magnitude of the clutter over the ocean is high it is recommended that data in the 
chirped mode be obtained over the forest due to the significantly lower backscatter 
associated with it at nadir. With these recommendations, at nadir, it is believed that 
rain rates > 5 mm/hr may be detected. The use of a better weighting function that 
gives lower sidelobe levels than that used (a raised cos 2 ) is also recommended. 

At off-vertical look angles all the range cells have a large clutter component 
associated with them due to the geometry. The use of higher angles of incidence (0 
> 60°) is recommended because of better SCR at these angles. With this 
recommendation, at off-vertical, it is believed that rain rates >10 mm/hr may be 
detected. 
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Various other techniques are described and recommended to improve the 
minimum detectable precipitation rate. These include trying to subtract the estimate 
of the clutter from the combined signal and clutter and trying to separate the Doppler 
of the rain echo and the surface echo. With these recommendations it is believed 
that it is possible to detect precipitation as low as 1 mm /hr at vertical and > 5 mm/hr 
at off-vertical look angles. 
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CHAPTER I 
INTRODUCTION 



1.1 INTRODUCTION 


The world ocean and water in and on land makes the climate of the earth 
uniquely suited for the sustenance of life. The process of evaporation and 
condensation of water vapor is essential for earth’s habitability. Rainfall is clearly 
central to water supply. Cloud radiation and fresh water from the rainfall may 
significantly affect the ocean circulation and marine biomass, while variations in low- 
level winds may determine up-welling and movement of nutrient-rich water (Simpson, 
1988). On an average, about 50% of the earth’s surface is covered by clouds but the 
ratio of the rain to the cloud area is seldom larger than 10% (Griffith et al., 1978). 
Although these estimates place an upper bound of 5% on the earth’s surface receiving 
rain at a given time, the actual percentage is much less (Barrett and Martin 1981). 
For the tropical zone between 30° N and 30° S, the occurrence of rain has been 
estimated to be about 4%, accounting for more than half the earth’s rainfall (Thiele, 
1987). Although contours of the global mean annual precipitation have been 
constructed over oceans and remote land areas, the annual rainfall at a particular 
location probably is known only within a factor of 2, at best (Thiele ,1987). 

One of the most sought after characteristics of precipitation is the rainfall rate, 
R,, defined as the volume of liquid water that falls through a unit area per unit time 
and normally expressed in units of mm per hour. Though rain events are dominated 
by lighter rainfall rates, higher rain rates account for most of the total liquid water 
reaching the surface. From approximate global distribution of rainfall, in the region 
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of highest rain rate, the of percentage time in a year the rainfall exceeds 50 mm/hr 
is 0.1%, or 526 min per year; it exceeds 150 mm/hr only for about 53 min annually 
(Crane and Blood, 1979). The contribution from the various rain rate categories to 
the total rainfall can be calculated from the probability density function (PDF) for the 
rain rate R,. The relative contributions of the various rain rate categories to the total 
rainfall are obtained by multiplying the PDF by the rain rate, R,, and integrating from 
zero to Rj. Available data show that the dynamic range required to characterize ocean 
precipitation in the tropics is less stringent than over areas with intense local 
convection. For example, in the GATE program, a radar that could accurately 
estimate rain rates between 2 mm/hr and 40 mm/hr would account for about 84 % of 
the total rainfall in a typical GATE area, but about only 64 % of the total rain volume 
over Texas, which has local convection (Chiu, 1988). However it appears that there 
is a correlation between the fractional area over which the rain rate exceeds a certain 
threshold on one hand and area-averaged rain rate on the other. If such relationships 
can be tuned to local climatology, for those applications in which area averaged rain 
rates are of primary importance, a large dynamic range may not be required. 

Most first-order weather service sites have the tipping bucket and the universal 
weighing-gauge precipitation monitors. The tipping-bucket gauges generally 
accumulate the number of 0.01 inch precipitation events in a day which they utilize 
to compare with the accumulation in other gauges. The rain rate estimate comes from 
the time between these tips. The universal weighing gauge is also capable of 
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providing rain rate information and is the main instrument utilized to provide the 5- 
minute to 1-hour precipitation data. By measuring the slope of the strip chart on the 
gauge, the rain rate may be obtained to at-least 5-minute resolution and even 1 -minute 
rain rates may be inferred from some charts. These charts are a good source of 
information for short-duration rainfall. Several investigators have developed more 
accurate techniques of estimating rain rate distributions directly from the rain gauges 
with resolution as low as 1 minute. 


1.2 MICROWAVE RAINFALL MEASUREMENT 

Rain drops, for many purposes, may be considered to be spherical under 
normal conditions (Battan,1973). Mie first devised a general theory that describes the 
scattering of a plane wave by a sphere. Hence, we are able to apply the Mie 
scattering theory to backscattering from rain drops with little approximation. For 
wavelengths under consideration (3 to 6 cm) it is usually possible to use the Rayleigh 
approximation (Battan, 1973). Various published results may be used to calculate the 
magnitude of the backscatter from the rain drops. Here, we are looking for an 
empirical relationship between the reflectivity / attenuation and the rain rate, to allow 
rain rate determination using any one of the available rain-retrieval algorithms. An 
accurate estimate of the dropsize distribution is crucial in determining these empirical 
relations. Observations of the dropsize distribution as a function of precipitation 
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intensity have been made in a variety of places and in various types of precipitation: 
widespread rain and snow from stratiform clouds, showery and orographic 
precipitation, etc. Because of the large variations observed in both time and space, 
the analyses of the data have been along statistical lines (Battan,1973). 

The dropsize distribution most widely used in the literature are the Laws- 
Parsons and Marshall-Palmer distributions, which compare reasonably well with one 
another. One of the major attributes of the Laws and Parsons dropsize distribution 
(Table 2.1) is that it covers a wide range of precipitation rates. Later chapters will 
discuss the theoretical relationship between the dropsize distribution and the empirical 
relationship and reduce it to a form that suits our purpose. 

Microwave rainfall retrieval techniques may be classified into 

1. Passive rainfall measurements 

2. Active rainfall measurements 

3. Hybrid rainfall measurements. 

1.21 PASSIVE MICROWAVE RAINFALL RETRIEVAL 

The ability of the microwave radiometer to extract rainfall depends heavily 
upon the contrast between the observed brightness temperatures over raining and the 
non-raining regions. The power emitted by a body in radiative equilibrium is a linear 
function of apparent or brightness temperature in the microwave region (Simpson, 
1988). 
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Due to the low emissivity (about 0.5) of the ocean, the ocean brightness 
temperature is very cold (< 150° K) (Ulaby, et al., 1986), and offers an ideal 
background to observe the emission from the rain at its equilibrium temperature. 
Land, in contrast has an emissivity of 0.7-0.9 which varies with vegetation, soil 
moisture, and other parameters. Thus, the brightness temperature measurements 
against warm and variable background appear to be inapplicable for retrieving 
atmospheric constituents such as rainfall. 

The absorption, emission, and scattering properties of the hydrometeors depend 
strongly on the frequency. Low frequencies are best described in terms of their 
absorption/emission effects, while high frequencies are best described by their 
scattering effects (Ulaby, et.al, 1981). Most recent passive rainfall retrieval 
techniques, however, rely on multichannel approaches that utilize the known 
scattering-to-emission ratios at the various frequencies to infer rainfall rate. 

The microwave brightness temperature depends on both the rain rate and many 
other atmospheric and cloud factors such as raindrop size distribution, cloud droplets, 
rain layer thickness, ice above the rain, inhomogeneities within the field of view, 
water vapor, wind speed at the ocean surface, and sea surface temperature. 
Multichannel techniques can retrieve some parameters along with the rainfall rate, but 
the behavior of the remaining parameters must satisfy some assumptions. 

Various algorithms exist that utilize the brightness temperature to estimate the 
rainfall rate, of these: 
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Decision-Tree algorithm (Kummerow, 1987) estimates the rain rate by examining the 
brightness temperature at several frequencies. Theory and observation confirm that 
the brightness temperature over the ocean tend to increase with the thermal emission 
for low rain rates, saturate and finally decrease due to the scattering from large drops 
and ice. A given brightness temperature at a given frequency will, in general, 
correspond to two rain rates. The examination of the brightness temperatures at 
several frequencies resolves this ambiguity. 

Dual-Polarization algorithm (Spencer, 1986) uses the polarization to resolve the 
ambiguity between low brightness temperatures produced by the precipitation and that 
produced by the ocean surface. The ocean surface, when viewed obliquely, produces 
radiance that is much more polarized than the precipitation. 

Multichannel Statistical methods (Olson, 1986) use statistical radiative-transfer models 
to extract the rainfall rate from the measured brightness temperature. These models 
allow the prediction of up-welling brightness temperatures for any combination of 
atmospheric and cloud parameters. Direct approaches minimize the difference in 
observed and calculated brightness temperature by varying the model. Statistical 
techniques use model-generated brightness-temperature/rainfall relations. 

1.22 ACTIVE-MICROWAVE RAINFALL RETRIEVAL 

Active rainfall retrieval algorithms are generally divided into three categories. 
At smaller rain rates, where attenuation is somewhat insignificant, algorithms based 
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on backscatter power work best. Attenuation based algorithms include the surface- 
reference technique (Meneghini, 1983), the dual-wavelength surface-reference 
technique (Meneghini, 1987), the standard dual-wavelength technique (Meneghini, 

1987) , the height-area rainfall-threshold technique (Atlas, 1988), and the mirror-image 
method (Meneghini, 1985). The least squares algorithm (North, 1988; Fujita, 1983; 
Kieu, 1992) estimates the rain rate profile by minimizing the difference in the profiles 
of the observed and estimated received powers. 

1.23 HYBRID MICROWAVE RAINFALL RETRIEVAL 

Attempts to combine radar and microwave radiometry retrieval have been 
successful. The calculated brightness temperatures from airborne radar data (Simpson, 

1988) agrees with the observed brightness temperatures. The multichannel linear 
rainfall technique described depends largely on the construction of cloud models that 
can accurately represent the observed brightness temperatures. With the correct 
model, results can be within 20 % for rainfall over an individual footprint (Weinman, 

1989) . 

The combination of the radar and the radiometer measurements, which places 
constraints on the radar equation, can yield a vertical hydrometeor profile from a 
single-wavelength radar. The two major constraints are the surface rainfall rate and 
the path-integrated extinction. The extinction coefficients near the surface can come 
roughly from the mean surface-rainfall rate obtained from the radiometer. The 
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extinction coefficients become boundary condition on the radar equation. The path- 
integrated extinction from the approximate hydrometeor profile derived from passive 
radiometry can be used as a constraint on the radar equation. 

Visible and IR radiometry is also widely used to confirm the presence of the 
clouds and measure their heights. The rain estimates in the VIS / IR radiometer are 
mainly from cloud observations. There is considerable interest in this technique 
because of the need for high spatial and temporal resolution of the precipitation 
events. Moreover, these methods can be used with geostationary satellites to obtain 
continuous coverage. 

Satellite IR radiometers can differentiate between the convective and stratiform 
rainy areas. The convective stratiform technique attempts to define the convective 
areas as those with relative minima in the IR temperature field, and attempts to screen 
out thin cirrus. The computed anvil background temperature determines the stratiform 
raining area. 

1.24 SPACEBORNE RAIN RETRIEVAL 

It is evident that the precipitation variability apparently coupled with the 
changes in the underlying surface ( particularly sea surface temperature) is associated 
with significant alteration in short-term climate, initiated in the tropics, which may 
lead to hemisphere-wide alteration in wind patterns and rainfall. Unfortunately the 
difficulty involved in accurately predicting this precipitation is an obstacle in 
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advancing understanding and prediction in this crucial portion of the earth system. 
Because of the inaccessibility of ground based sensors and the high cost of the surface 
radar networks, spacebome measurements offer the only means to obtain global 
rainfall data. 

The need for global precipitation data doesn’t in itself warrant the use of 
spaceborne radars. To justify its use requires that it be technologically feasible and 
cost effective in the sense that it cannot be duplicated by a more reliable or less 
expensive sensor. The unique features of the radar are well recognized. Unlike the 
lidar, it can penetrate through rain and cloud; unlike the radiometer it can vertically 
profile the rain and its sensitivity is not degraded by the high emissivity of a land 
background. While these benefits have long been recognized, in the absence of the 
necessary technology they are merely potential advantages. Without the advances that 
have occurred in the efficiency and reliability of the power amplifiers in the 
microwave and millimeter wave region, in low noise receivers, and antenna 
technology, spaceborne radars would suffer from poor sensitivity and reliability and 
limited scanning capability. 

One other reason for renewed interest in spacebome radars is the change in the 
attitude which views the radar as one among a set of complementary sensors on the 
spacecraft. The sensors on the orbiter, in turn, are treated as one of the many types 
of observation platforms, which include airborne radars and radiometers, ground based 
radars and rain gauge networks, and geosynchronous weather satellites. Although this 
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change in thinking doesn’t solve the inherent limitations of the spacebome radar it 
does provide a well defined role where its capabilities may be exploited with modest 
implementation. 

Although there is abundant literature on the measurement of precipitation from 
space using visible and IR sensors, microwave radiometry and radar ( Atlas and 
Thiele, 1981), none of it has been demonstrated adequately. Atlas and Moore (1987) 
suggest the use of spacebome SAR for rain measurements. Atlas et al. (1977) 
reported actual detection of precipitation by an airborne SAR. 

1.3 CHARACTERISTICS OF PRECIPITATION 

Rainfall measurements from spacebome platforms will improve our 
understanding of the hydrological cycle, atmospheric circulation, climatological 
models, and mesoscale rainstorm characteristics. However, no single model equation 
or set of parameters can be used to describe all forms of precipitation. To better 
understand and accurately apply the available models it is necessary to discuss the 
various precipitation types. For our interest we will be looking into the characteristics 
of: 

1. Stratiform precipitation. 

2. Convective precipitation. 

3. Hurricane precipitation. 
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1.31 CHARACTERISTICS OF STRATIFORM PRECIPITATION 

The life cycle of mesoscale convective systems consists of four stages: 

1. Formative, 2. Intensifying, 3. Mature, and, 4. Dissipating (Kieu, 1991). The 
convective cells dominate the total rain amounts in the formative stage (typically 2-6 
hr). During the intensifying stage, the stratiform precipitation begins to account for 
more of the total rain. In the mature stage, the stratiform component is equal to or 
greater than the convective component, as the total rain gradually decreases. The 
stratiform rain accounts for roughly 40% of the total rain during the lifetime of the 
weather systems. 

In the mid-latitude regions, stratiform rainfall typically has horizontal extent 
of hundreds of kilometers, durations exceeding one hour, and rain rates less than 
about 25 mm/hr. The spatial distribution of the total rainfall from one of these storms 
is fairly uniform. The rain rate averaged over several hours is usually similar for 
ground sites located up to tens of kilometers apart. This rain type usually occurs 
during spring and fall months and results, because of the cooler temperatures, in 
vertical heights of 4 to 6 km. 

In the stratiform region, the mesoscale updraft (Gamache, 1982) has been 
found between 20 to 40 cm/s above the altitude of the 0° C isotherm and a down-draft 
of nearly equal magnitude below this altitude. Using VAD (velocity-azimuth display) 
in the different portions of the stratiform area, the horizontal boundary between the 
mesoscale updraft and down-draft was found to occur just below the altitude of the 
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0° C isotherm (3.8 km). The magnitude of these mesoscale vertical motions in the 
stratiform region was 30 to 40 cm/s above the 0° C isotherm, and -10 to -30 cm/s 
below it (Chong, 1987). 

132 CHARACTERISTICS OF CONVECTIVE PRECIPITATION 

Convective rains arise because of the vertical atmospheric motions resulting 
in vertical transport and mixing. Strong vertical motions characterize convective cells, 
while weaker mesoscale motions are present in the stratiform region. The convective 
flow occurs in a cell whose horizontal extent is several kilometers across. The cell 
usually extends to heights greater than the average freezing layer at a given location 
because of the convective up-welling. The cell may be isolated or embedded in a 
thunderstorm region associated with a passing front. Because of the motion of the 
front and the sliding motion of the cell along the front, this high rate usually lasts 
only for several minutes. These rains are the most common source of high rain rates 
in the U.S and Canada (Kieu,1991). 

In studies of convective systems, maximum reflectivities exceed 45 dBZ*. The 
updraft maximum typically occurred at a height of 2.1 km. The average maximum 
updrafts are as high as 17 m/s. Other studies of convective systems (Chong, 1987) 
show that the maximum convective updrafts of up to 10 m/s were at 2.5 km along the 

Z is the reflectivity from the rain drops in the unit volume mm 6 /m 3 . 
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leading edge of the system, with a secondary maximum of 6-10 m/s above 5 km 
altitude 20-30 km behind the leading edge. The maximum convective down-drafts 
were 4 m/s. 

Radars operating at non-attenuating frequencies have been used to study both 
the horizontal and the vertical components of the convective rain systems. Crane and 
Blood report that rain rate variations of 100:1 could be observed over a range of 10 
km for a shower containing four intense cells. Similar measurements were made by 
Goldhirsh (1976), at Wallops Islands, VA. Goldhirsh (1976) also observed that the 
intense rain cells are elongated along the direction of motion. This direction 
correlated well with the average or median wind direction. 

1.33 CHARACTERISTICS OF HURRICANE PRECIPITATION 

Tropical cyclonic storms (Hurricanes) sometimes pass over the eastern 
seaboard during the August-October time period. These circular storms are typically 
50 to 200 km in diameter, move at 10-20 kilometers per hour, have extended melting- 
layer heights up to 8 km and have high ( greater than 25 mm/hr) rain rates. 

Experimental studies show the importance of predominantly stratiform 
precipitation surrounding the eyewall to the overall rainfall pattern of hurricanes. The 
eyewalls of hurricanes contain the heaviest rain rates, but encompass a relative small 
area. The precipitation in the surrounding region (out to 1 10 km), while having lower 
mean rain rates owing to its stratiform nature, covers much more area than the 
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eyewall. Up to 60 % of the rainfall within the inner region is predominantly 
stratiform precipitation in the region surrounding the eyewall. The eyewall rain rates 
increase as the storms intensify and the eyewall contracts. The increase in the mean 
eyewall rain rate is roughly balanced by the reduction in the eyewall area as the 
radius of the maximum reflectivity contracts. 

1.4 SIR-C / X-SAR PRECIPITATION EXPERIMENT 

The SIR-C/X-SAR precipitation experiment was proposed to study 
precipitation detection and measurement using the C- band channel of the Shuttle 
Imaging Radar (SIR) along with the X-band X-SAR. The prime objective of the 
experiment was to demonstrate that the spacebome (SAR) rain measurement concept, 
which would be extremely useful for the design and development of future global rain 
mapping missions. Although for climatological purposes rainfall must be observed 
over a long period of time, the results obtained from this proof-of-concept experiment 
will be useful to the studies of short-term rain phenomena and for verification of the 
radar-based rain-retrieval algorithms. Two other important objectives of the 
experiment were to establish the feasibility of measuring the rainfall motion from the 
radar Doppler signals, and to determine the modification of the ocean surface- 
scattering properties by the impinging rain-drops. 

It was expected that, with the SIR-C/X-SAR capabilities, quantitative rain 
backscatter measurements could be obtained in moderate to intense rain rate events. 
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By rotating the shuttle so that the antenna points to the vertical, it is possible to obtain 
the mean backscatter intensities from the rain-drops and the Doppler shift caused by 
the rainfall motion. With this observation geometry, it was expected to collect data 
that could demonstrate the use of the Synthetic Aperture Radar (SAR) technique to 
improve the along-track resolution. The SAR images obtained over the ocean during 
a rain storm could be used to study the changes in the surface characteristics due to 
the impinging rain drops. These results are expected to have important applications 
on the use of the ocean surface backscatter measurements for near surface wind 
estimation and for rain rate estimation using path-averaged attenuation. The rain data 
obtained by the SIR-C/X-SAR will constitute the first used and planned rainfall data 
set ever obtained from a spaceborne radar. 

The anticipated results of this experiment are 

• Quantitative evaluation of the rainfall rate. 

• Quantitative measure of the average rainfall velocity. 

• Quantitative measure of the height of the rain cloud. 

• Determination of achievable spatial resolution for rainfall measurement applications 

using synthetic aperture synthesis. 

• Quantitative measure of the modification in the ocean scattering mechanism as a 
result of the impinging raindrops, and comparison of such results with theory and 
previously obtained experimental data. 

• Polarization signature of the ocean and the signature’s variation in the presence of 
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precipitation. 


1.5 ORGANIZATION OF THE THESIS 

This thesis is based on the proposed SIR-C / X-SAR rain measurement 
experiment. After the experiment was accepted there arose a need for a feasibility 
study. This thesis presents the study and recommendations made for detection of 
precipitation from the spacebome platform. The thesis is organized as follows; 

Chapter 2 describes the concept of a synthetic-aperture radar and its ability to obtain 
better resolution compared to a real-aperture radar when the targets in the beam 
volume are correlated. However, when the targets have a random component the SAR 
fails in obtaining this fine resolution. It is shown that the resolution obtained by the 
SAR when the targets are not correlated is the same as that of a real aperture radar 
with the same parameters as that of the synthetic-aperture radar. 

Chapter 3 describes the various rain and clutter models. It also describes the 
relationship between the empirical rain models and the radar measured parameters. 
The rain rates are related to the radar reflectivity and attenuation by the power laws. 
The presence of the melting layer and an absence of sufficient knowledge of its 
electrical characteristics further complicates matters. However reasonable 
approximations can be made to account for it. Various ocean and land clutter models 


18 



at both X- and C-bands are considered, but only those of interest to us are presented 
here. The equations describing these models have been verified by comparison to the 
SEAS AT and ERS-1 data and with aircraft and land/tower-based data. These models 
are very close to experimental results. 

Chapter 4 discusses the geometry of the off-vertical viewing angles and calculations 
of SNR and SCR at these angles. The study was based on the SIR-C / X-SAR 
experiment and it was seen that there are two possible approaches to the rain rate 
retrieval. In the off-vertical mode the geometry is complicated because of the large 
beamwidth in the vertical direction. It is seen that all the rain cells are associated 
with a clutter component. This necessitates the calculation of an important 
benchmark, signal-to-clutter ratio (SCR). Though the SCR is important, in that it gives 
the minimum detectable precipitation rate, the Signal-to-noise-ratio (SNR) is also 
important. This chapter has both SCR and SNR calculations at off-vertical look 
angles. It shows that for off-vertical look angles the minimum detectable rain rate 
over the ocean is very high. To detect moderate rain rates ( >20 mm/hr) the angle 
of incidence must be high (> 60°). 

Chapter 5 discusses the geometry at vertical viewing angles and the SNR and SCR 
at these angles. Here at vertical it is assumed that the surface returns are more or less 
unaffected by wind speed, direction etc. Though the sidelobes of the antenna do not 
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affect near-vertical precipitation measurement, the sidelobes of the pulse do affect it. 
With the raised cos 2 weighting used, the range sidelobes are high enough to affect the 
SCR over the ocean adversely. However, over the forest, due to the low o° associated 
with it at vertical, there is a good possibility of detecting low rain rates ( > 5 mm/hr). 

Chapter 6 proposes that, theoretically at-least, it may be possible to detect low rain 
rates even in the presence of a large clutter component. This is an extension of the 
estimation of a low power signal in the presence of noise if there are enough 
independent samples of the signal and noise. It is shown that the number of samples 
required for successful estimation is within reasonable bounds. 

Chapter 7 discusses the assumptions made throughout the study and the conditions 
under which these assumptions may not hold. Also there are recommendations and 
comments about processing the rain data to accommodate these conditions. 
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CHAPTER H 

SYNTHETIC APERTURE RADAR RESOLUTION 



2.1 INTRODUCTION TO SYNTHETIC APERTURE RADAR (SAR) 

Synthetic aperture is a radar technique, implemented from an aircraft or an 
spacecraft, in which the effective synthesized antenna aperture is much larger than the 
real antenna. The radar motion is a prerequisite of the SAR technique. There are 
several ways to consider a SAR, each which lead to a different implementation of the 
processor. The synthetic aperture is so long that it is necessary to focus the beam to 
the point where the target is located, much as a camera is being focussed to near 
objects. In an analogy with the camera we may say that the lens focussed at infinity 
is the normal condition for an antenna. A lens focussed on a nearby object is 
analogous to the focussed synthetic antenna commonly used in the SAR. 

SAR is normally used to image the surface. Any moving object generate 
Doppler shifts that can compete with the Doppler generated by the radar motion; since 
the Doppler generated by radar motion is the key to SAR operation, it is difficult to 
overcome the problem of motion of the target in the beam volume. 

SAR signal processing is coherent phase compensation that corrects for the 
phase-shift history of the signal reflected from a specific point on the ground during 
the radar motion. With this compensation, only returns from that point are 
constructively added during the entire observation time. 
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2.2 SYNTHESIZED APERTURE 


To illustrate the concept of the synthetic-aperture, consider Fig. 2.1(a). The 
figure shows how one can build a real-aperture array of five elements focussed on the 
target T. Each element in the array is connected to a summing point, and all the 
voltages from different elements are added together at that point. The length of the 
transmission lines connecting each of the elements of the array to the summing point 
are based on the distance to the target and the position of the particular element in the 
array. The phase delay for each signal arriving to the summing point from the 
element must be the same if the total contribution is to be added in phase. When the 
signals add in phase the array is said to be focussed at that target point. This means 
that the equivalent total distance from the target to the summing point must be the 
same for every element, so that 

Li + R) — 1^2 + R2 = L3 + R3 = L 4 + R 4 = L s + R5. 

The lengths Lj are the "equivalent free-space lengths " of the transmission line (i.e. 
lengths if the phase velocity in the transmission line is the same as that for the wave 
in space). During transmission, the signals all originate from the summing point for 
reception, so that during transmission the phase shift along each path is the same as 
that along every other path. Consequently, for the location of the target T, the total 
round-trip phase shift is the same for each element, and the electric fields from the 
different elements all add in phase. For any other point in front of the antenna, the 
distances between the target and the summing point differ, so the contributions from 
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the different array elements differ in phase, and the output is not as much as it is for 
a target located at the focal point T. In actual real-aperture antennas, the lengths L, 
to L 5 may not be as indicated but, may differ by a wavelength or several wavelengths 
in such a way that the phase relationship is maintained. 

We may extend the concept of the real-aperture equivalent focussed array for 
a moving antenna by storing the signal voltages from the target T in memory 
locations as shown in Fig. 2.1(b). Here we assume that both the amplitude and the 
phase of the arriving signal are recorded. After all the memory locations are 
occupied, the contents of these locations are phase shifted by amounts needed to 
compensate for the differing distances R, through R 5 . In essence we have multiplied 
the voltage V, by e 2jkR i. The phase correction 2kRj accounts for the round-trip phase 
shift travelled by the signal from element i. The phasor sum of the signals is 5 times 
the value for one signal. If the same phase correction had been applied to signals 
received from different points, such a phase shift in space would be different than the 
correction, and the magnitude of the resultant phasor would be less. Hence, applying 
proper correction results in focussing the synthetic array at the location T. 

In this implementation of the SAR there is no requirement for continuous 
motion from array element 1 to 2 to 3 to 4 to 5. The elements need not be uniformly 
spaced, and the motion may be such that the different element positions are not even 
occupied in any particular order. The only requirement is that the position of each 
element relative to the focal point T be known and the appropriate phase correction 
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applied in the processor. 

Mathematically we can state this description of the total phase shift on path 
1 of Fig. 2.1(a) as <|> T1) the phase shift in space on that path as <t) Sh and the phase 
shift on the transmission line as <(> L ,. Thus we may write 

4>n = 4>li + <1 >S2 = 4>ri = 4> L i + <1 >5i = 4>r (2.1) 


Here the phase shift in space is 


4>si ■ 2 k R, 


( 2 . 2 ) 


and the phase shift on the transmission line is 

= 2 k L L i 


(2.3) 


where k L is the wave number for the transmission line. 

The required phase shift in the transmission line is 

4>Li =4 > r " 4>si 


(2.4) 


From (2.4) and (2.2) the length of the transmission line can be determined as 

r _ <t> T - 2 KR ± (2.5) 

1 ' 

The voltage received at the element i is therefore just 
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( 2 . 6 ) 


V = e ~ j2kSi = e ~ j * si 
Therefore, we need a correction 

4>C2 = 4>Li = 4>r “ 4>si ( 2 . 7 ) 

If we apply the correction, the resulting voltage at the summing point is 

V = ~ j * r = 5 e j * T ( 2 . 8 ) 

i 

(for the target at the focal point). For the synthetic-aperture it is usual to make 
<(> T = 0 by making <() Ci = - <J> Si . 

Various schemes have been devised to perform true synthetic-aperture 
processing. In these processes one accumulates all the signals associated with 
different elements of the synthetic array and simultaneously multiplies with a 
correction factor. The processing may be performed in real time by multiplying each 
return at a particular range by different correction factors. The results of the 
multiplication can be stored in different locations corresponding to different synthetic- 
apertures being constructed. 

To relate the synthesized aperture to the resolution, consider the geometry in 
the Fig. 2.2. The radar antenna on the airborne vehicle is shown in three positions 
A,B and C. At location A the forward edge of the beam first intercepts the target; at 
position B the beam is abreast the target and at position C the rear end of the beam 
is leaving the target. Thus the total possible length of the synthetic-aperture L p is the 
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(a) Antenna Beams Illuminating Target at T 



(b) Two Adjacent Elements in Synthetic Array 


Fig. 2.2 General idea of the geometry of the Synthetic aperture 
(Ulaby ct al„ 1982) 
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resolution that would have been obtained by the real antenna in a real-aperture SLAR. 


Lp=V h rRo < 2 * 9 > 

Here (3^ is the horizontal or along track beamwidth of the real antenna, the subscript 
r is used to distinguish this beamwidth from the beamwidth of the synthetic-aperture 
antenna. 

In part b of Fig. 2.2, two adjacent elements in the synthetic array separated by 
a distance Ax, are shown, with portions of the ray R a and R b going from these 
elements to the target. Since the elements are close together, the rays may be 
considered essentially parallel. Hence, the difference in the distance along the two 
rays is 

R a - R b =Ax sinP (2.10) 

For the real aperture antenna the phase shift associated with this distance is 

A <t> r = K( R a - R b ) = k Ax sinP (2.11) 

For the synthetic-aperture array, the phase shift both as the wave goes from antenna 
to the target and as it returns from the target to antenna must be considered. Thus in 
this case 

A<t> s =2K{R a - R b ) =2 k Ax sinp (2.12) 

For an antenna of length L ( real or synthetic), where L is long enough that the 
beamwidth is small, the value of p associated with the first null is given by solving. 
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for the synthetic-aperture. 


2kL sinP - 2kL p = ti (2.13) 

and for the real aperture, 

kL sinP - kL P = tc (2.14) 

Hence the value of P associated with the first null for the synthetic-aperture is P s = 
n / 2 Kl, whereas for the real aperture it is p r = n / Kl. This factor of two difference 
occurs at all angles associated with the pattern for the synthetic-aperture. Hence the 
beamwidth for the synthetic- aperture is half that of the real aperture having the same 
length. 

The along track beamwidth P hs for the synthetic-aperture is 

ft = (2.15) 

Whs 2 L 

where a^, is the aperture illumination taper factor associated with the synthetic- 
aperture. 

Using the expression for the beamwidth of the synthetic-aperture, the along- 
track resolution is obtained as 

r = Q h R = — a hc (2.16) 

This expression applies regardless of the actual length of the synthetic-aperture used. 
If the entire potential length L p is used, the finest possible resolution r ap for the SAR 
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is given by 


_ ^ Ra hs _ _ a /is ( 2 . 17 ) 

3P ‘ 2 L p 2 (J hr 

Thus, the resolution attainable is inversely proportional to the beamwidth of the real 
antenna used in the SAR. 


Piir ~ a hr 


X 

1 


( 2 . 18 ) 


where a hr is the aperture illumination taper factor for the real antenna and 1 is its 
length. This value of p hr may be used in (2.17) to obtain the potential synthetic- 
resolution in terms of the length of real antenna and the aperture taper factor as 

_ * 1 a hs „ J_ ( 2 . 19 ) 

ap 2 a hl X 2 

when the illumination taper factors are the same. 

Thus the finest possible resolution is just half the length of the real aperture. 
It is independent of the range and also the wavelength. The resolution r ap is the finest 
possible resolution that one can attain with the synthetic-aperture radar having an 
antenna fixed to a spacecraft that is flying in a straight line, and for a stationary 
target. 
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2.3 RESOLUTION OF THE SAR 


The range resolution achieved by the SAR is the same as that of the 
conventional side looking real-aperture radar (SLAR). The range resolution for the 
SAR is 

r D = C (2.20) 

R 2 

where c is the speed of light and x p is the is the rectangular pulse width. The 
resolution on the ground is 

r = c t p (2.21) 

y 2 sin6 

where 0 is the angle of incidence. 

The resolution possible by the SAR in the azimuth direction is more complex. 
To obtain an expression for the azimuth resolution, let us consider the SAR as a 
focussed Doppler-beam- sharpening processor (see Fig.2.3). In Fig 2. 3. (a) the radar 
is on the aircraft at a distance x R away from the origin and at a height h. The 
coordinates of the radar are (x R , 0, h). The coordinates of the target are (0, y T , 0). 
The ellipse shown represents the half power contour on the ground of the beam of the 
real antenna carried by the aircraft. The solid outline surrounding the target and 
bounded in the x-direction by the edge of the beam, shows the area contributing to 
the return at a particular instant of time. Also shown are the contours of constant 
Doppler frequencies (isodops) surrounding the target, with spacing Af D . 
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A bandwidth B Df of a Doppler filter is used for beam sharpening (Wiley, 
1967). The Doppler frequency for the target in this case is given by 

-2 u {x R - x T ) 

XR ( } 

The target location x T is fixed, whereas the radar location x R is varying. For the 
illustration the Doppler frequency is positive because x R is negative. The Doppler 
frequency decreases linearly with the passage of the beam over the target, if we 
neglect the slight variation of the range. 

Parts (b) - (f) of the figure illustrate the effect of the focussed Doppler beam- 
sharpening system as the target is successively illuminated first by the forward edge 
of the beam and later by other portions of the beam, until finally the beam passes it. 
The along-track length of the illuminated area on the ground is simply the resolution 
that would exist for the real-aperture radar having the same beamwidth. 

r ar = P b R (2.23) 

The total bandwidth of the Doppler-shifted signal that is received by the radar is the 
same as would be found for a real-aperture radar signal. It is given by 

= X P * = 'Tl 5 ' <2 ' 24) 

where X is the wavelength of operation, R is the slant range to the target, u is the 
velocity of the aircraft P h is the horizontal beamwidth and r ar is the azimuth resolution 
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of the real-aperture. 

The Doppler beam-sharpening is accomplished by using a narrow bandwidth 
filter that tracks the Doppler-shifted frequency from the target. The limits on the 
ground defined by the filter are shown by the hatched region of width r a . Its 
bandwidth is indicated as B Df . 

The Doppler frequency shift obtained by the beam passing by the target (see 
parts (b) - (f) of Fig. 2.3) is shown in part (g) of the figure. The Fig. 2.3(g) shows 
the limits of the tracking filter passband following the Doppler shift associated with 
the point target. The Doppler frequency from the other parts of the beam are 
discriminated against because the Doppler frequencies associated with them cannot 
pass through the tracking filter. The bandwidth of the filter controls the along track 
resolution possible. From (2.24) the possible along track resolution is 

r = A* B d£ (2.25) 

a 2 u Df 


2.4 SAR RESOLUTION IN RAIN 

SAR can use the full synthetic-aperture only when the targets in the pulse 
volume are correlated, i.e., their phase positions remain essentially fixed. This is 
because the coherent phase compensation corrects for specific phase shifts from a 
point, only returns from that point are constructively added. To find the maximum 
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Fig. 2.3 Doppler beam sharpening SAR. Concept assumes a narrow-band 
Doppler filter. (Ulaby et al., 1982) 
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possible aperture that can be synthesized during T a , the time during which the target 
is illuminated, let us consider the bandwidth between the isodops for a SAR. 


*df = 


2 

R A 


( 2 . 26 ) 


where U is the velocity of the platform, R is the range and the others have their usual 
meaning. The ratio between the synthetic-aperture length L s and the time T a to pass 
that length is simply the platform speed. Thus the maximum possible aperture that 
can be synthesized during the integration time T a is 

L s = U T a ( 2 . 27 ) 

A larger Doppler spread (a vs ) leads to shorter integration time. For a smaller 
synthetic-aperture where the resolution element is large. Atlas and Moore (1987) show 
that the Doppler bandwidth is 


B df {precip) = 4 


( 2 . 28 ) 


From (2.28) we may write 


r as ^ 


„ R 

1 / 


( 2 . 29 ) 


and also 


35 



L=k — — (2.30) 

4 o 

v vs 

and the synthetic-beamwidth is 

B = 2 (2.31) 

5 R 

All this means that the synthetic-beamwidth varies and is broadened in proportion to 
the Doppler spectral width of the target. Further, it means that the number of pulses 
that can be coherently integrated and the effective length of the synthetic-aperture L s 
are inversely proportional to a vs . 

This applies only when 

°vs <-£ L Z n <2 - 32) 

where 

U = the velocity of the spacecraft in m/s, 

L n = length of the along-track real aperture, 
a vs = standard deviation of the Doppler velocity spectrum, 

For X-SAR the standard deviation of the Doppler velocity (due to the wind 
turbulence, shear, etc.) must be less than 3.75 m/s for a three-fold improvement in the 
resolution for the SAR, compared to the RAR (i.e., for the synthetic-beamwidth to be 
one-third of the real beamwidth). 

This indicates that the effective beamwidth increases with increases in a vs , the 
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spread of the target Doppler spectrum. Further, the integration gain decreases with 
a vs because the number of pulses that can be coherently integrated along track is 
proportional to the target coherency time and thus inversely proportional to ct vs . 
Therefore, the echo power is independent of the resolution. The echo power is 
identical to the conventional radar with an along-track beamwidth equal to the real 
beam used in the SAR. The SNR equation for the SAR looking at the rain will be 
the same as that for the real-aperture radar except when the SAR prf exceeds that 
required for Nyquist sampling of the Doppler spread. (Atlas and Moore, 1987) 
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CHAPTER III 

MODELS USED IN THE CALCULATIONS 



3.1 INTRODUCTION 


To determine the signal-to-noise ratio (SNR) of the signal returned from rain, 
the scattering coefficient of the rain must be known. Though theoretically it is 
possible to model the rain as discrete scatterers (apply Mie scattering theory with 
Rayleigh approximation), it is not possible to describe all the different rain types by 
a single equation. This is in part due to the drop distribution of the volume scatterers. 
Here, in this study, for the SNR and SCR calculations to be of any use rain rate must 
be the principal variable. It is easier to use empirical relationship developed by 
meteorologists for these calculations. These empirical relationships can be linked to 
the Rayleigh approximations by use of model dropsize distributions. 

The SCR calculations require knowledge of the surface scattering coefficient 
in addition to the scattering coefficient due to rain. The rain-measurement experiment 
is to be primarily over ocean and forest. Although various mechanisms can be used 
to describe the ocean backscatter, here we must select an ocean-surface echo model 
that shows variations with wind speed, angle of incidence and aspect angle (direction 
of wind with respect to the radar boresight). 

Tropical forest clutter at vertical is unrelated to wind, and varies little with 
time of day or season. The actual measured scattering coefficient over the Amazon 
at Ku band rather than a model, will be used in the SCR calculations (Ulaby, et.al, 
1986) (Schroder et al., 1985). 
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3.2 BACKSCATTER FROM THE OCEAN 


Radar measurements of backscatter from the ocean were first made under 
controlled conditions during World War II. At that time, and for many years 
thereafter, the measurements were made to increase knowledge about the nature of the 
sea clutter that obscures radar echoes from surface targets and submarine periscopes. 
The NASA earth-resource program in 1964 gave impetus to the radar observation of 
the sea for research. This was followed by a series of airborne experiments leading 
to the use of radar scatterometer to measure windspeeds and culminated in the Skylab 
S-193 and Seasat space-based scatterometer experiments. 

Radar backscatter from the ocean at most angles of incidence is strongly 
related to the speed and direction of the wind. When the wind starts to blow over the 
calm ocean, the first waves formed are due to frictional drag. As these short waves 
build up, non-linear interactions transfer energy to waves with larger amplitudes and 
longer wavelengths. This interaction and transfer continues until the equilibrium is 
reached. For angles of incidence above 20° the predominant radar scattering 
mechanism is Bragg scattering, while at near-vertical incidence angles the mechanism 
is a quasi-specular one in which geometric or physical optic techniques can be used. 
The amplitude of the radar signal beyond 20° is dominated by the Bragg-resonant 
ripples and the local angle of incidence. 

The variation of the scattering coefficient (o°) with wind speed for a particular 
radar look direction is reasonably well described by 
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o° = A W 


(3.1) 


where U is windspeed in meters per second and A and y depends on the polarization, 
angle of incidence (0), look direction, and frequency. Expressing (ct°) in dB, (3.1) 
can be written as 

a°{dB) = 10 log o° = 10 log A + 10y log ( U ) (3.2) 

An orthogonal regression procedure is used to fit a straight line relation between (o°) 
and log (u). Values for the exponent y are obtained for different radar and ocean 
parameters. 

3.21 SURFACE ECHO MODEL 

To obtain the surface-echo model we decided to use data from the wind-vector 
measurements. The use of the radar scatterometer to determine the surface wind 
vector over the ocean is a proven technique. The Skylab S-193 Radscat instrument 
allowed measurements over much of the world ocean surface, but there were very few 
concurrent surface observations. The measurements were compared with an 
underflight by the AAFE Radscat instrument, and the values of the a 0 observed by the 
two instruments were within 0.5 dB of each other. This indicated that the calibration 
of the two instruments were very close and were accurate in the absolute sense 
because they were calibrated independently. Fig. 3.1 and Fig. 3.2 are data from this 
experiment. 
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Fig. 3.2 Regression fit for Skylab a 0 measurements at near vertical incidence. 

Note the difference in response near vertical and away from vertical 
shown in Fig. 3.1. 
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The AAFE Radscat was the best calibrated airborne instrument used prior to 
1980 in measurement of radar sea returns. Experiments with this instrument were 
aimed at improving the knowledge of the wind speed and directional response. The 
data were analyzed by both the NASA Langley Research Center group and workers 
at The University of Kansas and the City University of New York. 

The circle-flight experiment, developed by Jones et al. (1978) at NASA 
Langley Research Center, allowed the accurate measurement of azimuth variation of 
echoes from the sea for the first time. An example of such measurement is shown in 
Fig. 3.3. Previous measurements were only for up, down, and cross winds. A basic 
model function for azimuth variation was developed from the Jones measurement: 

o° = A + B cos<|> + C cos2<|> (3.3) 

A comparison with theoretical results of Chan and Fung (1977) is given in Fig. 3.4. 
Theory and experiment correspond, allowing for the variation in measurements. 
Because of the excellent calibration of the AAFE Radscat, and because there were 
many flights accompanied by high quality surface observations, this was one of the 
most reliable aircraft wind-response data sets available until recently. 

Soofi (1978) developed an empirical model to calculate the magnitude of the 
ocean returns. This model matches closely with the AAFE Radscat, SEASAT and 
SKYLAB data (Schroeder, et al., 1985). It also accounts for various wind speeds and 
aspect angles (direction of the wind with respect to the radar boresight). 
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Fig. 3.3 Example of the circle-flight observation of the azimuthal variation of 
a 0 at 13.9 GHz, along with the regression fit (Moore, et al., 1978) 
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Fig. 3.4 Comparison of AAFE Radscat 13.9 GHz a°measurement with theory 
of Chan and Fung (1977)(Moore, et al., 1978) 
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The equation for this ocean-clutter model is 


o° db (f,0, W, <|>) = A + Be os <|> + C cos 2<(> (3.4) 

where the coefficients A, B and C take values depending on the polarization used. 
For the polarization of interest (VV) the values of A, B and C are functions of wind 
speed, angle of incidence and frequency. 

A = (-24.642 + 10.8256 logW ) + (-0.1983 - 0.0017 logW ) 0 
+ (1.17 25 - 0.6191 logW ) f + (-0.0123 + 0.0124logW r ) f 0 

B = (0.4337 - 0.7041 log*/ ) +(-0.0085 + 0.0143 logW ) 0 
+ (-0.0082 + 0.0153 logW ) f + (0.0004 - 0.006 log AT ) f0 

C = (0.1619 + 1.4160 logW ) + (-0.0086 + 0.0089 logW ) 0 
+ (0.0067 - 0.0124 log2 ) f 

Here the frequency is in GHz, the wind speed in m/s and the angle in degrees. This 
model predicts the value of the coefficients within about 2 dB. The maximum error 
occurs at 4 GHz. 

Schroeder, et.al (1982) proposed the SASS (SEASAT-A Satellite 
Scatterometer) model based on data at 13.9 GHz. This is an empirical function giving 
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the relationship between the ocean normalized radar cross section and the wind vector 
at a height of 19.5 meters above the surface. The relationship is defined in the form 
of a table (Table 3.1(a)) which gives the two coefficients G and H in the equation 

NCRS(dB) = 10 [G(0, 4>) + H(6,4>) log 10 U] (3.5a) 

where 0 is the radar incidence angle, <)) is the angle between the wind direction and 
the radar look direction, and U is the neutral stability wind speed in meters per 
second. The variation of the NRCS (c°) with the surface wind speed obtained from 
Soofi’s and the SASS model are compared in Fig. 3.5. Though the model proposed 
by Soofi is good for preliminary calculations, it was never in the form of published 
results. The models are close for upwind direction and downwind direction, but the 
crosswind values differ. 

For C-Band there were recent models available from ESA based on the ERS-1 
experiment. ESA’s CMOD4 (Ad Stoffelen, et al., 1992) gives the model as a 
complex transfer function. The transfer function is given by 

o 0 = K b [l +lc x cos(<j>)+Jc 3 tanh ( k 2 ) cos (2<|>) ] 16 I0 (a+Y ♦ <V+ P >) 

where a, P, y are coefficients in Table 3.1(b) (Ad Stoffelen, et al.,1992 ) and 
Legendre polynomials of the first and second order and <j> is the direction with respect 
to the beam, k,, k 2 , k 3 are functions of the coefficients and the wind velocity v. The 
Legendre polynomials have 0, the angle of incidence as their variable. 
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Schroeder Soofi 


Fig. 3.5 Comparison of the scattering coefficient (a) of the ocean surface 

obtained from the SASS (Schroeder) and Soofi’s model at 9.7 GHz, 
60°, cross wind. 
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Preliminary calculations of the backscatter at C-band based on the CMOD4 
model show an unusually high y. Soofi’s model of backscattering from the ocean is 
inconsistent with this. However, Soofi’s C-band model is based on interpolation 
rather than actual experiment. Further, CMOD4 is in the form of published results, 
based on the ERS-1 experiment. Hence the CMOD4 model is preferred over Soofi’s 
C-Band model; a comparison between the two is shown in Fig. 3.6. 

3.22 EFFECT OF RAIN ON SURFACE RETURN 

The modification of the backscatter due to impinging rain drops may be 
important. In the area of the ocean where it is raining, the cross section increases 
due to the impinging rain drops. Moore, et al. (1979) showed that, for a given wind 
speed, the backscatter cross section depends on the rain rate, as shown in Fig. 3.7. 
Hansen (1986) later conducted both laboratory and field experiments that show a 
significant enhancement of the echo intensity as a result of the splashes and ripples 
due to the impinging raindrops. This enhancement of the radar return applies only for 
low wind speeds. At high wind speeds the effect of rain rate on backscatter 
modification is lost; i.e., the backscatter at these high wind speeds is independent of 
the rain rate (Moore, et al., 1979; Blivens and Norcross, 1988) as seen in Fig. 3.7. 
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Fig 3.6 Comparison of the scattering coefficient (a) of the ocean surface 

obtained from the CMOD4 and Soofi’s models at 9.7 GHz, 60°, Cross 
wind. 
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Fig. 3.7 Experimental results indicating the loss of backscatter modification of 
the ocean surface due to impinging rain drops, at higher surface wind 
speeds. (Moore el al., 1979) 
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3.3 BACKSCATTER FROM THE RAIN 


Mie first devised a general theory that describes the scattering of a plane wave 
by a sphere. Mie theory shows the backscattering cross section of a spherical drop 
(Battan,1973): 


°b = 


7 ia‘ 


S (-D* (2/77+1) ( a B -b m ) 

/n = l 


( 3 . 6 ) 


where a is the drop radius and <> equals 2rca /A. The coefficient a m and b m are in terms 
of spherical Bessel functions and Hankel functions of the second kind, with arguments 
X and n c , the complex refraction; 


n c = n - in 


( 3 . 7 ) 


where n is the ordinary refractive index and K is the absorption coefficient of the 
material involved. For spheres that satisfy the Rayleigh condition ( % «1), the 
backscattering cross section of a single particle is proportional A. -4 . That is, (Battan, 
1973) 

a b = — ! K\ 2 D 6 ( 3 . 8 ) 

b A 4 

where K is related to the complex permittivity n c and is approximately 0.92 when the 
scatterer is liquid water, over the wavelength band most commonly used in 
meteorology (Table 3.2). Summing over N v particles per unit volume, the rain 
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volume backscattering coefficient a v is expressed as 


° v = S 0*0^) 

i-l 

-5 "v 

= JL_ SD 6 (3.9) 

* 4 i - 1 

We can write Zj D 6 = Z, the reflectivity, which is usually given in mmV. 

Radar meteorologists have derived an empirical power law, the Z-R 
relationship, from the available data. That is, 

Z = a, r :/** (3.10) 

where a, and b t have been found to depend on the dropsize distribution, rain type, 
rain rate, frequency, and region. The coefficients differ markedly at small and large 
values of R r . However, with the exception of orographic rains, most differ only 
slightly at rainfall intensities between 20 and 200 mm/hr. Although it is clear that 
the Z-R relation is useful, the values of the coefficients a, and b { have to be selected 
judiciously. A single Z-R relation for all types of rain cannot account for the 
variations for different types of rain (Table 3.3) Battan (1973) calculated the Z-R 
relation using the Mie scattering theory. Table 3.4 shows the various relationships for 
different wavelengths and dropsize distributions. Clearly for a given wavelength, but 
for different rain rates, a, takes drastically changing values (Table 3.4). However, for 
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modest and small rain rates, assuming the value of a, to be around 300 and b { to be 
about 1.5 is reasonable. Here we use the relation 


Z = 300 R^ 5 


(3.11) 


to describe the rain reflectivity. The scattering coefficient of the rain is described by 


o y - * 5 300 \ k \ 2 

v A . 4 


(3.12) 


when X is in centimeters and rain rate R in mm/hr. 


3.31 ATTENUATION IN THE RAIN VOLUME 

The attenuation experienced by the radar returns from the rain may be 
attributed to two effects: 1. Absorption, and 2. Scattering of the power away from the 
radar beam (Ulaby, et.al, 1981). This is because, in the cloud or rain volume, the 
scatterers are usually assumed to be randomly distributed within the volume. Thus, 
there are no coherent phase relationships between the fields scattered by the individual 
particles. 

Mie scattering theory leads to two cross sections; the total scattering cross- 
section Q s and the absorbing cross-section Q a . The extinction cross-section Q, is the 
sum of Q s and Qa- The extinction cross-section is the area which, when multiplied 
by the incident intensity, gives the total power removed from the plane wave by the 
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particle. The absorption cross section is the area which, when multiplied by the 
incident intensity, gives the power dissipated as internal heat in the particle. 
Atmospheric gases, cloud droplets, and Rayleigh raindrops act as absorbers, but both 
scattering and absorption are notable for non-Rayleigh (large) raindrops. 

Assuming that the Rayleigh condition is valid (i.e.,% « 1), the Mie scattering 
and absorption cross sections reduce to the corresponding Rayleigh scattering and 
absorption cross sections. That is (Ulaby, et.al, 1981) 

r 2 

k s = fp(r)Q s (r) dr ( 3 . 13 ) 

where k s is the volume scattering or absorption coefficient, p(r) is the dropsize 
distribution, Q(r) is the scattering or absorption cross-section of the sphere of radius 
r in m 2 , r, and r 2 are the upper and lower limits of drop radii contained in the cloud. 
Expressing (3.13) in terms of scattering efficiency = Q/Jtr 2 and the dimensionless 
parameter % = 2kt/X the scattering coefficient gives 

x 2 

k s = -——r f X 2 P(X> 5 S (%) dx ( 3 . 14 ) 

8rc 3 J 

Xi 

where the limits on the integral cover the entire range of possible values of X- 
Expressing (3.13) in terms of absorption efficiency ^ = QJnr 2 and the dimensionless 
parameter %, the absorption coefficient is 
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(3.15) 


= “rr/x 2 P(x) 5 a (x) <*x 

8 71 3 J 
Xi 

For a given dropsize distribution the volume extinction coefficient for (= 

K + K) is 


x 2 

k er = P(%) ^ er (x) cfX (3.16) 


Here, as previously defined, % = 2ixr/X. and ^ er (X) * s the Mie extinction coefficient for 
rain given by 

i ei (n lX ) = 4 2 (2/n+l) i?e[a m +2?J (3.17) 

X 2 »-i 

Here Re signifies the real part of and b m , the Mie coefficients (Ulaby, et.al, 1981, 
PP 291). 

From the practical point of view it is desirable to relate the extinction 
coefficient for rain k er directly to the rain rate R r . The relationship commonly used 
in literature (Gunn and East, 1954) is of the form 

= K, R? (3.18) 

where k, is in dB km' 1 per mm hr' 1 , Rf is the precipitation rate in mm hr' 1 and b is a 
dimensionless parameter. The coefficients K, and b, which are functions of the 
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wavelength X, are determined by fitting (3.18) to computed values of k„ for the 
known dropsize distributions. The dropsize distributions in turn determine the 
precipitation rate through the relation 

R r = 6ti 1CT 4 I,v i d\ (3.19) 

i«i 

where v^ms' 1 ) is the terminal velocity and dj(mm) the diameter of the i* rain drop and 
N v (m' 3 ) is the total numbers of drops per unit volume. Table 3.5 (Ulaby, et.al, 1981) 
provides list of values for K, and b at each of the several microwave frequencies. 


3.4 RADAR EQUATION FOR METEOROLOGY 

The basic radar equation for any target is 


P t G 2 A. 2 a 
(4n) 3 i? 4 


(3.20) 


where the variables have their usual meaning. 

If Z Cj is the average total backscatter cross section of particles per unit 
volume (m 3 ), the radar cross section may be expressed as 


o = V m Eo i (3.21) 

where V m is the volume of the resolution cell. We can use the Rayleigh 
approximation when the dropsize is « X to write. 
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(3.22) 


2 




Tt 5 2P 6 

A. 4 


K 


■ 2 
I 


where is A. in centimeters and K is a factor involving the complex permittivity. Using 
(3.12) we may write the radar equation as 


= P c G 2 Tt 2 V m 3oQ zl 5 , k]2 (3.23) 

64 A, 2 i? 4 

This radar equation does not include the various system losses, nor does it account for 
attenuation due to clouds, rain, the ice melting layer and the atmosphere. Additional 
factors added to (3.23) can account for this. 


Pr = 


P t G 2 


* 2 Ye 11 h 

(4tt ) 3 R 4 


Pi 


10 -o.i* #r 


(3.24) 


is the final form of the radar equation for meteorology obtained by substituting 

additional factors to account for losses. 

where 


n 


it 5 300 i? 1 ~ 5 I JTl 2 

X 4 


IQ ' 10 


(3.25) 


G is the gain, B is the bandwidth, L s is the system loss, L^, is the signal loss due to 
absorption by the ice melting layer, t| is the rain reflectivity, (3.18) is the 
attenuation due to rain. The other parameters have their usual meaning. 

Dividing (3.24) with the noise power we can calculate the SNR. The volume 
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of the rain resolution depends upon the angle of incidence. These volume 
calculations, based on the geometry at different angles of incidence, are given in the 
next two chapters. 
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CHAPTER IV 

CALCULATIONS AT OFF- VERTICAL ANGLES 





Surface, for 1 pulse length 


Fig. 4. 1 (b) Off-vertical look geometry for the plane earth case used to calculate the 
pulse volume. 
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4.1 INTRODUCTION 


For vertical incidence and a narrow transmitted pulse, the backscattered power from 
the hydrometeors arrives at the receiver prior to the large surface contribution. At 
non-vertical viewing angles, however, the surface return arrives simultaneously with 
the desired echo (see Fig. 4.1(a)). Thus, the surface clutter can mask the returns due 
to the rain at the surface. To calculate the SNR and the clutter coupled to the 
resolution element, the volume of the scattering element and the corresponding surface 
area coupled to it must be calculated. For the SAR the volume encompassed by the 
radar beam is similar to that of a real aperture radar. Thus each resolution element 
in the radar beam will be similar to that of the real aperture radar. 

When the radar is pointing off-vertical, some of the resolution elements in 
the rain volume have different characteristics. This is in part because of the finite 
height of the rain volume. It is possible to identify different resolution cells in the 
rain volume based on whether they are completely or partially rain filled. 

4.2 GEOMETRY 

The geometry shown in Fig. 4.1(b) is for non-vertical viewing angles making 
the plane-earth assumption. Here, for simplicity in approximate calculations, we 
assume a sharp beam edge and neglect the sides of the mainlobe and side lobes. More 
exact calculations would require accounting for true beam shape, but these simplifying 
assumptions enable us to define range cells that are sharply beam limited. Thus we 
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may divide the beam into three well defined regions and allows computation of the 
dimensions of the observed volume as a function of the look angle. This assumption 
may not hold true in reality but SCR calculations are made on cells that are close to 
ideal. The beam below the cloud top is divided into range cells corresponding to a 
slant-range resolution of 15 m. The swath R s on the ground for vertical beamwidth 
(3 V and satellite height h is 

R s = h [ tan (0,) - tan(0)] ( 4 . 1 ) 

where 

e, = p v + e (4.2) 

and 

R S1 = (h - hj [tan (0,) -tan(0)] ( 4 . 3 ) 

is the swath width at the cloud top, and h, is the thickness of the rain. Computation 
of the swath gives us an idea of the coverage possible by the X-SAR mission. 

The number of rain-filled cells is 


where 


n 


_ -^lmax 


( 4 . 4 ) 


•^lmax 


cos (0 1 ) 


( 4 . 5 ) 
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off nadir geometry illustrating the three regions ol intrest 


Geometry at off-vertical look angles, with completely rain filled 
volume, identifying three regions. 






is the range in the rain-filled volume and r R is the range resolution. Fig. 4.2 shows 
three cases of cell lengths. 

Region 1. Pulse position at the end of the return that extends up to the edge of the 
beam, not to the rain top. 

Region 2. Pulse position in the center of the return that extends from the ground to 
the rain top. 

Region 3. Pulse position at the start of the return that extend from the inner edge of 
the beam to the cloud top and includes no ground return. 

If we assume that the rain extends throughout the swath, most of the returns 
correspond to those in region 2. Further, since we can calculate the signal-to-noise 
ratio and the signal-to-clutter ratio using the real-aperture parameters, the width of a 
cell in the direction normal to the plane of the paper, for power calculations is 

r An - P a [R - nr R ] ( 4 . 6 ) 

where P a is the azimuth beamwidth, R is the slant range to the ground, r R is the range 
resolution and n is the cell number, measured from maximum R backward. 

4.21 CELL VOLUME CALCULATIONS 

In region 1 (longest range) the volume of the n* cell from the bottom is 

V n =Lr R r a ( 4 . 7 ) 

where 
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L = n r R cot (0 1 ) 


(4.8) 


is the length of the rain cell in the plane of incidence and normal to the ray from the 
radar to the rain. The width of the rain cell in the direction normal to the plane of 
the paper is r a , r R is the slant-range resolution and 0 av is the average angle of 
incidence. 

In region 2 (center of the beam) the volume is the same for all pulses. 




V - r r 

n a * sin (0. J 


(4.9) 


Here 0 av , the average angle of incidence, is used when calculating the range to the 
pulse volume (bounded by concentric spheres) to avoid the complication arising 
because it varies slightly across the beam. The ground area coupled to the rain cell 
in the above two cases is 


A - r t 


sin(0 r ) j 


(4.10) 


where 0 r is 0 t in region 1 and 0 av in region 2. 

In region 3 (shortest range) there is no ground return. The length of the cell 


is 


L = m R r cot (0) (4.11) 

where m is the index starting from the top. The cells in region 3 are valid only for 
the assumed beam and minimal antenna sidelobes. 
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Even if the rain cell fails to fill the swath completely three similar regions may 
be identified, but the number of constant volume rain cells is less, as seen from 
Fig. 4.3. 

Region 1 The region for which returns come from the far edge of the rain cell. 

Region 2 The region for which returns come from near the center of the cell. 

Region 3 The region for which returns come from near the edge of the cell. 

The volume of the rain cells at the center is the same as that in region 2 of the 

completely filled beam. However, for region 1 , if 


h n = n 


c °s(0 air ) y 


and the length of the cell is 


L = 


sin (0 a „) 


and the volume of the cell can be written as 


( 4 . 12 ) 


( 4 . 13 ) 


V n = r R r a 


n 


cos(0 av ) sin (0 air ) 


( 4 . 14 ) 


where n is an integer and h^ is defined as the distance to the pulse intersection with 
the outer edge of the cell. 

In region 3 


h m =m £r — - 

l cos (0 av ) 


( 4 . 15 ) 
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where m is an integer value corresponding to the cell number (starting from the top) 
and h m is the height to the intersection of the range cell with the inner edge of the rain 
cell. The length of the cell from this point to the top is 


L = 


sin (0 av ) 


( 4 . 16 ) 


The volume of the rain cell in this region is 


v m = r a r R 


( h - 

{ sin (0 av ) 


which reduces to 


( 4 . 17 ) 


V m " ® r a 


rl 


sin(0 av ) cos(0 av ) j 


( 4 . 18 ) 


Thus the volume of the resolution cell in the rain is a function of the look angle and 
the range resolution only, since r a is assumed to be a constant. 

For plane-earth geometry, the look angle and the incidence angle are the same. 
However, for the spherical Earth they are not. Since the rain height is small 
compared to the satellite height, we can use incidence angle in the calculations (i.e., 
we may use all the equations discussed previously with incidence angles rather than 
the pointing angle). 
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4.3 SIGNAL-TO-CLUTTER RATIO 


The backscatter method uses the returned echo from the rain to compute the 
rain rate. For non-vertical angles of incidence, however, the surface returns arrive 
simultaneously with the desired echo and increase in magnitude as the range gates 
within the mainbeam begin to intersect the surface (see Fig. 4.1(a)). Thus, the surface 
clutter can mask the return due to the rain. For the backscatter method to work, the 
echo from the rain should be substantially above the surface echo to allow separating 
them. Hence, we compute the signal-to-clutter ratio (SCR). This is defined as the 
ratio of the desired signal (rain return) to the clutter signal (surface return). Since 
the surface and rain are excited by the power from the same transmitter, one need not 
use the radar equation in its entirety for these calculations. The scattering coefficient 
for the rain has the units m 3 but the scattering coefficient of the ground has the units 
m 2 . To compare the returns, one must multiply the rain returns by the volume of the 
rain cell and the surface returns by the area of the surface coupled to the 
corresponding rain cell. 

The scattering coefficient for the rain from (3.22) when X is in cm is (an 
additional factor 10 10 accounts for change in the units) 

O v = -Hi Z \K\ 2 10- 10 (4.19) 

A 


The volume of the rain cell is given by 
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(4.20) 


V 


r R r a 


L 


where r R is the range resolution, r a , the azimuth resolution and L is the length of the 
rain cell, discussed in the previous section. Since the scattering coefficient of the 
ground calculated from the model is in dB, we can get the total cross section by 
adding 10 log(ground area) to c dB , obtaining a gdB as 

a gdB =°dB + 10 log (r a R sl ) (4.21) 

where r a and R s] are from geometry discussed in the previous section. 

To determine how much the rain signal lies above the ground signal we compute 

Ao <* - o vlda - o ffdB ( 4 . 22 ) 

where 

°vj d8 = 10 log (o r ) (4.23) 

and 

o r = o v V (4.24) 

Now let us introduce the concept of ideal rain cells. An ideal rain cell is 
defined as a resolution cell that has a completely rain-filled volume. Clearly these are 
the kind of rain cells that are of interest for our calculations. From the previous 
discussion it is clear that the rain cells in region 2 come close to our definition of an 
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ideal rain cell. 


4.31 SIGNAL-TO-CLUTTER RATIO IN THE MAIN PART OF THE RAIN 

Let us now consider a more realistic case, including attenuation effects. For 
non-vertical look angles, if the entire beam volume is assumed to be rain filled, all 
the returns from the ground are attenuated. Further, as we move from the top to 
bottom of a rain cell, the signal is attenuated increasingly as shown in Fig. 4.4. To 
consider this type of attenuation for each rain cell (only using Region 2), we must 
determine the attenuation factor for different path lengths. Thus the actual surface 
echo is 


°s = 


o c e" B 2 Rlnax 


(4.25) 


where R lmax is the range in the rain. The actual rain echo is given by 

°r = r R r a L o v k e[ (4.26) 

where 


k 


GZ 



f -a dR 

e * 


dl 


(4.27) 


This accounts for the different ranges for the same rain cell. 
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Fig. 4.4 Increasing attenuation of the returns from the same rain cell as we 
move from the top to the bottom. This is due to the longer path 
traversed through the rain volume by the returns from the bottom 
compared to those from the top. 
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The rain return involves 


Rmax 

- j a dR (4.28) 

o v e * V m 

v m 

where V m is the volume and c v is the scattering coefficient and a is the attenuation 
coefficient in nepers/km. Since the range changes along the same rain cell (see 
Fig. 4.4) it is necessary to integrate the attenuation over these ranges. We can write 
(4.28) as 


Vgff X a Z R 


-« / dR 

o v | e * du 


(4.29) 


From Fig. 4.4 


R m Kax ~ 

^max ^nax 


(4.30) 


which may be written as 


R = 


_ [ *^max ^ 


^nax 


R max 


(4.31) 


Now substituting = h, / cos6 and u max = h, / sin0 in (4.30) we have 


r = (i- u sin 8 ) 




COS0 


(4.32) 


and substituting in (4.32) in (4.29) we have 
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a ^tnax 

°v v eff = r a r R o v e " 55 * 5 f e“ u tan ® du (4.33) 

0 


Evaluating the integral, and making the substitution u roax tanG = h, / cosG we have 


V 


eff 


r a 

a tanG 


(l - e 


a hi 
cos0 j 


(4.34) 


Now, if we write 

V eff = r a r R L e ff and compare both sides, we see that 


L 


eff 


1 - e- a *™ 
a tanG 


(4.35) 


The equation is used to account for the variation of the attenuation in the same rain 
cell. This works only for a large enough that L is not limited by the beamwidth. 
If L be the actual length of the cell we may write L eff = K^L , where is the 
attenuation factor. Evaluating 


( 1 _ o ® ^ ^lmax \ 

k er = * COS(0) (4.36) 

which is the two-way attenuation of each rain cell. Here R lmax is the range in the rain 
volume, and oc from Table 3.5, is the attenuation coefficient in nepers/km, using the 
dB-to-nepers conversion. Substituting (4.36) in (4.26) will give the attenuated rain 
echo; similarly subtracting (2 a dB R lmax ) from (4.25) will give the attenuated surface 
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Fig. 4.5(a) 


t r 

10 

Rain rate in mm / hr 


100 


-X- 6 m/s 24 m/s 


Plot of the SCR in dB vs Rain rate in mm/hr at X-band, 60°, cross 
wind, for the given wind speeds. 
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Fig. 4.5(j i>) Plot of the SCR in dB vs Rain rate in mm/hr at C-band, 60°, cross 

wind, for the given wind speeds. 
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SNR in dB 



Fig. 4.6 Plot of the SNR in dB vs Rain rate in mm/hr at X and C bands for look 
angles of 30° and 60°. 
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Tt 5 300 i? 1 ■ 5 ! Jc ! 2 
A. 4 


( 4 . 39 ) 


*1 


IQ-10 


G is the gain, B is the bandwidth, L s is the system loss, L,, is the signal loss due to 
absorption by the ice melting layer, r\ is the rain reflectivity, k er is the attenuation due 
to rain from (4.28), and F is the receiver noise figure. The other parameters have 
their usual meaning. The plots in the Fig. 4.6 are made from the data generated using 
(4.30) for different angles of incidence at X-and C-band. 

This SNR calculation is based on the assumption that the entire beam volume 
below the cloud top is rain-filled. If this is not true the SNR will decrease. The 
result of this is that the minimum detectable rain rate increases for fixed SNR. 
Studies show that a typical rain storm diameter for rain rates of 10 to 15 mm is 10 
to 15 km (see Fig. 4.7) (Goldhirsh and Rowland, 1982). Thus, relative large storm 
diameters may be expected at lower rain rates, and the probability of partial beam 
filling along track is low. If however the cross-track beam is only partially filled we 
may consider the region 2 discussed earlier to a be special case of regions 1 and 3; 
i.e., there will be at least a few ideal rain cells. 
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Fig. 4.7 Typical rainstorm diameter for the given rain rates (Goldhirsh and 
Rowland, 1982) 
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CHAPTER V 

CALCULATIONS AT NADIR 



5.1 INTRODUCTION 


In Chapter-IV it was concluded that it would be possible to measure rain rates in 
excess of 30 mm/hr at angles of incidence greater than 60°. Here we consider the 
radar pointed at vertical. At vertical it is seen that, due to the non ideal nature of the 
transmitted pulse, there is considerable clutter coupled with the rain echo. When 
measurements are made at vertical, the problems are somewhat different than at the 
angles well away from vertical. Further, the geometry is different because of the 
presence of both beam-limited and range limited conditions. This means that the 
volume of the resolution cell depends upon its position in the beam. 

5.2 GEOMETRY 

The echo power for the SAR is identical to that for the conventional real 
aperture radar, using the same antenna. Hence, will use the same geometry as that 
of real-aperture for power calculations. Fig. 5.1 shows the geometry for the vertical 
incidence. The rain volume may be divided into three regions: 

Region 1. where there is partial beam filling due to the pulse intercepting the cloud 
top. 

Region 2. where the cells are assumed to be completely rain filled. 

Region 3. where the rain cells intercept the surface. 

Here we consider calculations only for those cells that are completely rain 
filled (ideal rain cells). 
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Fig. 5.1 Geometry at vertical illustrating the volume of one pulse and 
identifying three regions. 
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B 


Fig. 5.2 Geometry at vertical on the 
of the surface coupled to ea 


8 . 



return, where a dB is the attenuation dB/km. We substitute these attenuated returns in 
(4.22) to obtain a realistic estimate of the SCR; this is plotted in Fig. 4.5(a) and 4.5(b) 
for X- and C-bands. 

It is seen that for low rain rates the X-band signal lies far below the clutter for 
lower angles of incidence of (0 < 50°). However, for higher angles, the clutter lies 
below the signal for moderate rain rates at low wind speeds. At C-band the clutter 
is very large. Fig. 4.5(b). 


4.4 THE SIGNAL-TO-NOISE RATIO 

From Chapter-Ill we have the radar equation for meteorology. From Fig. 4.2 the 
volume of the radar beam for the off-vertical look angle may be expressed in terms 
of the horizontal and the vertical beam widths for a pulse of duration T as 

v m jx Ph R hi c t (4.37) 

m 4 2 sin (0 av .) 

where h, is the thickness of the rain region. This volume corresponds to the cells of 
interest; i.e., the ideal rain cells discussed in the previous sections. The resulting 
equation for the SNR is 


SNR = 


ft G 2 X 2 c x hi T) L s Lfc 
K T B F 512 n 2 i? 3 sin(0 av ) 


(4.38) 


where 
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Fig. 5.3 



Geometry to calculate the total number of range cells cluttered with the 
surface echo. 
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The volume of these rain cells is 


v n = -J P* P v r j? i?n 2 (5.1) 

where R n is the range to the cell in question, (3 h and p v are the horizontal and the 
vertical beam widths, and r R is the vertical resolution. The factor of n/4 accounts for 
the elliptical shape of the beam. The ground area coupled to each cell except the first 
can be calculated, assuming the area is approximately rectangular (see Fig. 5.2), as 

A n * 2 Ph R n ( ~ x(n-l) ) (5.2) 

and the ground area coupled to the first cell is given by 

A i = 2 Ph R n x(D (5.3) 

where 

x(n) = R n sin (0 n ) (5.4) 


6 n = cos' 1 ( — ^ ) 

K n 


(5.5) 


and 


R n = h s + n r R (5.6) 

We define h s as the shuttle altitude, R n as the range to the particular cell in question 
and n as the cell number ranging from 1 to m. To calculate m, the total number of 
cells coupled to the ground, we use the relation 
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m = 


R s - h s 


(5.7) 


where R s is the maximum slant range to the ground (see Fig. 5.3) 


Rs cos ( p v /2 ) 


(5.8) 


5.3 SIGNAL-TO-NOISE RATIO 

From Chapter-Ill, the radar equation for meteorology for vertical look angle 
is 


SNR n = — 

n K T B F 


G 2 X 2 p 


h c t 1) L b L s k 


1024 7t 2 R 2 


[Rr] 


(5.9) 


where 




7T 5 300 i? 1 - 5 lie! IQ- 10 

X i 


(5.10) 


Here G is the gain, B is the bandwidth, L s is the system loss, 1^, is the signal loss due 
to absorption by the melting ice layer, r| is the rain reflectivity, F is the receiver noise 
figure, and the other parameters have their usual meaning. 
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From Fig. 5.4 it is seen that the volume occupied by the rain above each rain 
cell changes from one end of the cell to the other. This is due to the fact that pulse 
volume is curved whereas the rain top is not. Hence the signal from the cell edge 
suffers less attenuation than the center. To account for this it is necessary to integrate 
the attenuation over the entire pulse volume. Here it is necessary to make the 
assumption that the cell has constant attenuation in the azimuth direction. The rain 
returns involve 

K 

2 

a v -1 r* f R 2 e' a R * dR (5.11) 

_K 

2 

we may write R = (h + n r R ) where n is the cell number under consideration, r R is the 
range resolution and h is the height to the shuttle from the rain top. To account for 
the variation of the distance travelled by the signal in the rain volume it is assumed 
that the beamwidth can be divided into number of small angles. From Fig. 5.4 

R *- h ' nI * "cos® (5 - 12 ' 

where 0 varies from - Pv/2 to Pv/2. Substituting (5.12) in (5.11) for the rain returns 
we have 



K 

2 ( h (co3fl - l) 

{h + n r R ) 2 j e cos ^ 

Pv 


+ n 


de 


(5.13) 
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Here the attenuation coefficient 




d 0 


(5.14) 


k er [Rr] 


JU 

2 

/■ 

.h 

2 


-2 a [fir] (* tf ( -£22^11) ♦ /.*) 


d0 


(5.15) 


where 


a [J?r] = b_ Rr^ (5.16) 

8.66 

The values of b and c depend on the frequency used (see Table 3.5), h d is the shuttle’s 
height above the top of the rain cell and h k is the distance of the rain cell from the 
rain top. Using (5.9) we plot the SNR in dB vs rain rate in Fig. 5.5 for X- and C- 
bands. 

Only calculations for the cells that are entirely rain filled are presented here. 
For partially rain-filled cells, the SNR decreases in proportion to the rain-filled 
volume. 
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SNR in dB 



Fig. 5.5 Plot of the SNR in dB vs Rain rate in mm/hr at X- and C-Band for the 
vertical geometry. 
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5.4 SIGNAL-TO-CLUTTER RATIO. 

The backscatter method uses the returned echo from the rain to compute the 
rain rate. For vertical incidence the surface returns, coupled by the antenna sidelobes, 
arrive simultaneously with the desired echo and increase in magnitude as we approach 
cells near the surface. Since the echo from the rain should be substantially above the 
surface echo, we compute the SCR. 

The scattering coefficient of the rain, when X is in centimeters, is 


a v = — Z \K\ Z 1CT 10 
A. 4 


(5.14) 


Using the appropriate volume from (5.1) and the extinction coefficient from (5.15) 


o r = o v V 


k gr [ Rr ] 

Pv 


(5.18) 


and the appropriate area coupled to the rain cell from (5.2) in 


= o 0 (A n ) , 


(5.19) 


where the magnitude of <f depends on the type of surface. We compute the SCR 


A o rfn = o r - o„ 

™ r dB 9dB 


(5.20) 


where 


° rdB = 10 log (o r ) (5.21) 

and 
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a gdB = 10 lo gio^) 


(5.22) 


For the ideal case, the clutter coupled to the rain returns by range sidelobes is 
negligible. However in practice the surface return from range sidelobes arrives at the 
receiver at the same instant as the desired rain echo. The range sidelobe level with 
a raised cos 2 window is low, with the first side lobe about 21 dB down. However, 
due to the very large magnitude of the surface clutter and due to the fact that the side 
lobes fall too slowly (see Fig. 5.6), almost all rain cells in the present configuration 
have significant clutter coupled to them. Fig. 5.6 shows that the magnitude of the 
sidelobes that are coupled to the surface echo depends on the cell for which the 
calculations are made. Further, the sidelobes that are coupled to the ground return 
change in magnitude with the cell under consideration. To account for this, using 
(5.7) and the cell number under consideration, it is possible to obtain the exact 
number of side lobes coupled to the ground. Their magnitude Sl dB (n), can be obtained 
from the figure. The weighted clutter coupled to the cell under consideration is 
obtained as 


<V <*r> = 2 o g ( R z ,n ) i 0 °- lsJ <® (n) (5.23) 

n 

where 


a g {R r ,n) 


A n K or [Rr ] 10 01 °° 

Pv 


(5.24) 


Since the vertical-incidence surface echo is so strong, returns due to the antenna 
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sidelobes are negligible compared to those due to the range sidelobes. Hence all 
clutter calculations are made using the range sidelobes. 

Near vertical incidence the signals produced by the Bragg-scattering 
mechanism are dominated by those produced by the geometric optics or physical- 
optics mechanism. In this region the slopes are such that one can achieve specular 
reflection from the facets on the surface of the sea, and many facets having narrow 
backscattering patterns are close enough to the specular condition that signal returned 
from them by the quasi-specular mechanism are strong. At and near vertical 
incidence the strongest reflection would occur with a perfectly flat surface or one that 
was very undulating. As the surface gets rougher, more of the vertically incident 
energy is scattered away from the radar receiver. Hence the radar backscattering 
decreases with the increasing wind-speed and wave height near vertical, whereas it 
increases with increasing wind speed at angles beyond 12°. 

In vegetation scattering, the air-vegetation boundary usually is unimportant. 
Hence, in general we have a combination of volume scattering in the vegetation layer 
and surface scattering by the underlying ground. It is seen that the surface scattering 
decreases rapidly with increasing 0 for angles close to vertical, while the volume 
scattering decreases very slowly with increasing 0. Unless the canopy is dense we 
expect that the backscatter due to the soil to be at-least as important as the volume 
scattering. In tropical rain forests, the canopy is dense enough that we need not 
consider the backscatter from the soil. 
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For vertical incidence the scattering coefficient of the ocean at 13.9 GHz is 
close to 12 dB for wind speed of 10 m/s (Schroeder, et al.,1985) and that of forest to 
be -6 dB (Ulaby, et al.,1986). Using (5.20), plots of the SCR for various rain 
thicknesses over ocean are compared with SCR over forest for rain thickness of 1 km 
at X-band (see Fig. 5.7(a)). 

The transmitted pulse is chirped and the weighting used is a raised cos 2 with 
a 0.4 pedestal. Ideal cos 2 weighting is difficult to achieve because of the imperfection 
in amplitude and phase response and various non-linearities in amplifiers. However, 
even if the ideal cos 2 weighting could be achieved it would widen the range 
resolution. Fig. 5.7(b) shows the comparison between the SCRs obtained using a 
raised cos 2 window and an ideal cos 2 window over the ocean at X band. Fig. 5.7(c) 
and 5.7(d) show corresponding comparisons at C band. From the figures it is seen 
that the X-band and C-band SCR with a raised cos 2 window is poor over the ocean 
for shallow rain thickness. However over the forest the SCR appears to be reasonable 
even for shallow rain thickness with a raised cos 2 window. Further for both X and 
C bands the SCR is expected to be good even for shallow rain thickness over ocean 
if an ideal cos 2 weighting could be achieved. The X-SAR has the capability of 
transmitting a short 3 psec pulse and processing it without dechirping it (Fuk Li, 
1994). This results in a wider range resolution (450 m) but the sidelobe problem is 
eliminated. This means that rain at an altitude as low as 500 m above the ground can 
be detected. 
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SCR in dB 



Rain Rate in mm 


over ocean over forest 


Fig. 5.7(a) Comparison of the SCR in dB vs Rain rate in mm/hr over forest to that 
over ocean at X-band for the vertical incidence. 
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SCR in dB 



raised cos~2 ideal cos A 2 


Fig. 5.7(b) Comparison of the SCR in dB vs Rain rate in mm/hr over the ocean 
at X-band, vertical incidence, obtained by using the raised cos 2 
weighting to that obtained by use of ideal cos 2 weighting. 
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SCR in dB 



Fig. 5.7(c) Comparison of the SCR in dB vs Rain rate in mm/hr over forest to that 
over the ocean at C-band, vertical incidence. 
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Fig. 5.7(d) 
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Rain Rate in mm 
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Comparison of the SCR in dB vs Rain rate in mm/hr over the ocean 
at C-band, vertical incidence, obtained by using the raised cos 2 
weighting to that obtained by use of ideal cos 2 weighting. 
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CHAPTER VI 

ESTIMATION OF THE SIGNAL IN THE 
PRESENCE OF LARGE CLUTTER 



6.1 INTRODUCTION 


The principal limitation for the use of a spacebome imaging SAR as a rain 
radar is the surface-clutter problem. Signals may be estimated in the presence of 
noise by averaging large numbers of independent samples. We applied this method 
here to obtain an estimate of the rain echo by averaging a set of N c samples of the 
clutter in a separate measurement and subtracting the clutter estimate from the 
combined estimate. 

Chapters IV and V show that there is adequate signal-to-clutter ratio (SCR) for 
the experiment to be feasible only at high rain rates that have a small probability of 
occurrence. However, signals from smaller rain rates that are more common are 
swamped by the surface clutter. The model used in the SCR calculations at X-band 
is the SASS model described by 

NCRS(dB) = 10 [G(0, x) + H(Q,x) log 10 U] (6.1a) 

and the model used in the SCR calculations at C-band is the ESAs CMOD4 described 
by 

o 0 = K b [1+iq cos (<j >) +k 3 tanh(ic 2 ) cos(24 >)] 1 - 6 io ♦ (V+ P> > 

Fig. 4.5(a) and (b), which have plots of SCR vs rain rate for off-vertical angles of 
incidence, show that, at lower rain rates, the SCR (in dB) is negative. This is also the 
case for vertical incidence, as seen from Fig. 5.7(a) and Fig. 5.7(c). This situation is 
analogous to the estimation of a low signal in a noisy environment. It is possible to 
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estimate the signal even when the SNR (in dB) is negative, if many independent 
samples are available. We extended this concept to estimating the signal with 
negative SCR. 

6.2 MEASUREMENT IN THE PRESENCE OF CLUTTER 

A radiometer operates with a signal-to-noise ratio (SNR) that is much less than 
unity, and radars with a reasonable number of independent samples of fading can 
operate with an SNR of the order of unity or even less. Two samples may be 
considered as independent if they are far enough apart so that the correlation between 
them is essentially zero, which usually occurs when the spacing between the samples 
is large compared with the reciprocal of the bandwidth. Each pixel in a side-looking 
airborne radar (SLAR) image usually represents an average of several independent 
samples, although the number is small. Each synthetic-aperture radar (SAR) pixel 
could represent only one sample, but it usually contains 3 or 4 samples (multilook). 

At the input of the receiver the received voltage is the sum of the mean values 
of the signal v/ and the noise v n . 

/ ( 6 . 2 ) 

= Vj + 

The mean value of the received power is the sum of the mean values of the signal 
and noise power because the random phase causes cancellation of the cross terms in 
the mean, giving 
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( 6 . 3 ) 


P = P ' + P 

r i r s c n • 

where P n is the noise power. 

Now, the signal power may be broken up into two components P/ = P s + P c , 
where P s is the power returned from the rain and P c is the power returned from the 
surface. We may express the mean received power in terms of the SNR S n and SCR 
S c as 


= Pjl 




( 6 . 4 ) 


When N r independent square-law-detected samples of the received signal are averaged, 
the standard deviation of the average is 


a rdN ~ 



( 6 . 5 ) 


When the SNR is high enough the estimate of clutter alone can be obtained by 
averaging a set of N c samples of clutter in a separate measurement. The standard 
deviation of the clutter alone is (neglecting noise) 


°cdN ~ 


Pc 

c 


( 6 . 6 ) 


An estimate of received signal can be obtained by taking the difference between 
estimates of combined signal and clutter and clutter alone. 


103 



( 6 . 7 ) 


~ ^rd Pcd 

The variance of either a sum or difference of two independent random variables is the 
sum of the variances of the individual components, hence the variance of signal alone 
is 


Os = ® rdN + ® cdN 


( 6 . 8 ) 


Now if we assume that compared to the SCR, the SNR is very high, and substitute 
from (6.5) and (6.6) we obtain , 


- S / , _L, C 

O s ~ — (1 + — ) + — — 


N, 




( 6 . 9 ) 


Rearranging 


(lt t )2 


f ± )2K 


( 6 . 10 ) 


N 


N r 


where K = N r / N c . When N r = N c K is 1. We may use (6.10) as a measure of the 
estimate of signal in presence of clutter. Solving for the number of independent 
samples N r , for K=l, we have 


N r 




2 



( 6 . 11 ) 


At lower SCR the number of samples required for successful estimation 
increases rapidly. 
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Fig. 6.1(b) Number of samples required to estimate the signal within 20% at 60°, 
X-Band for 6 m/s (60X6) and 24 m/s (60X24) cross wind. 
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Fig. 6.1(c) 
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number of samples 



ocean forest 


Fig. 6.2 


Number of samples required to estimate the signal within 20% at 
vertical, X-Band. 
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For example, (6.1 1) shows that a change in SCR of 3 dB from -10 dB to -13 dB gives 
a fourfold increase in the number of samples required for successful signal estimation. 
Fig. 6.1(a) and 6.1(b) give us a bound of the required number of samples for different 
ocean-surface wind speeds at angle of incidence 60° at X band. Fig. 6.1(c) has the 
required number of samples for successful estimation of the signal for a wind speed 
of 6 m/s at C band. These plots give us an idea of the number of independent 
samples required to estimate the signal to within 10 % to 

20 % accuracy at off-vertical angles of incidence. Here, for rain rates less than 10 
mm/hr, the number of samples required is very large. This is because the SCR at 
these rain rates decreases to small values and the numerator of (6.11) increases 
rapidly. Fig. 6.2 has plots of the number of independent samples required at the 
vertical angle of incidence for different rain heights over the oceans and 6.2(b) has 
the same parameters at C band. As in the off-vertical case the number of samples 
required for successful estimation for rain rates > 5 mm/hr is very large, but feasible. 

6.3 POSSIBLE NUMBER OF INDEPENDENT SAMPLES 

Two approaches may be used to obtain independent samples in a SAR. In one 
approach the finest possible resolution r ap is obtained and the resulting fine resolution 
image has these pixels averaged together. In another approach several looks are taken 
from slightly different directions, each with the final resolution r aa , and these pixels 
are averaged together. In terms of the observation time we may write the number of 
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independent samples as 


N r = A f D T 


2 Ph R 


( 6 . 12 ) 


which reduces to the simple relation 


N r 

r ap 


( 6 . 13 ) 


6.31 OFF- VERTICAL ANGLES OF INCIDENCE 

For the X-SAR the dimensions of the antenna are 12.1 x 0.3 m. The 
potential resolution is 6.05 meters. We write the actual azimuth resolution as 


r 


aa 


2 O v R 

U 


( 6 . 14 ) 


where (3 h is assumed to be r aa / R. Values of a v may go up to 6 m/s. The actual 
resolution, for an altitude of 215 km, r aa is 688 meters at 60°, assuming 6 m/s 
turbulence. The number of samples available in azimuth is about 113. If we obtain 
N r independent samples in the range by degrading the range resolution, the total 
number of independent samples available is 113N R . The effect of degradation of 
resolution on the SCR is minimal because any change in the surface area is offset by 
a corresponding change in the rain volume. 

From Fig. 6.1(b), for the angle of incidence 60° and rain rate 10 mm/hr, the 
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number of samples required to estimate the clutter within 20 % is 30000 for a surface 
wind speed of 24 m/s. We have 1 13 samples in the azimuth and need about 300 in 
range. This number in range can be obtained by degrading the range resolution to 
4500 m from the possible 15 m. For lower wind speeds the required number of 
independent samples is much smaller. 

At 60°, C-band, for wind speed of 6 m/s (see Fig. 6.1(c)), using similar 
calculations, it is seen that the number of samples required to estimate the signal is 
enormously high. Hence the minimum detectable rain rate is still > 10 mm/hr. 


6.32 VERTICAL INCIDENCE 

The potential resolution of the X-SAR is 6.05 meters. We write actual the 
azimuth resolution as 


r 




2 a v R 

U 


(6.15) 


where p h is assumed to be r aa / R. Values of o v may go up to 6 m/s. The actual 
resolution, for an altitude of 215 km, r aa is 344 meters at vertical. The number of 
samples available in azimuth is about 57. By degrading the range resolution and 
hence the range sidelobes, the total number of independent samples available is 57 N R . 
However, since the number of rain cells available are limited in this case compared 
to off-vertical angles of incidence, the range resolution cannot be degraded to an 
arbitrarily large value. 


Ill 



( 


Fig. 6.2 shows that, for a rain rate of 5 mm/hr, estimation within 20 % and a 
rain cell 3 km high, the number of samples required is 323. We have 57 samples in 
the azimuth and need about 6 in range. This number in range can be easily obtained 
by degrading the range resolution to 90 m from the possible 15 m. For a range cell 
2 km high the corresponding numbers are 83 1 samples required with 1 5 in range and 
a degraded range resolution of 225 m. For a range cell 1 km high, however, the 
required number of samples is very large, and it is not possible to obtain a reasonable 
estimate for resolution cells below this height. 

6.4 ESTIMATION OF THE SIGNAL 

To obtain an estimate of the signal using this method, first a measurement of 
the signal and clutter is made from the rain volume. Then, with the same resolution, 
measurement of the clutter is made outside the rain volume. To obtain an adequate 
number of independent samples, the degradation of resolution may be performed 
during post-processing by averaging N R x N A samples. 

The accuracy of the estimation heavily depends on the knowledge of the radar 
backscatter cross section in the presence and the absence of the rain. The ratio of the 
surface backscatter radar cross section in the presence of rain <f to that in the absence 
of rain (f , is assumed to be unity. This, however, is a questionable assumption since 
the cross section is modified due to the impinging rain drops (Hansen, 1986). To 
overcome this limitation it is suggested that the estimation of the clutter be made in 
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the same rain volume but with the synthetic aperture. When the synthetic aperture is 
used, because of the narrow Doppler bandwidth, the interference due to rain is 
minimal. However, the rain appears as a noise component to the clutter estimate. 
This noisy clutter estimate is assumed to give better results than that obtained using 
the clutter estimate made outside the rain volume. 
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CHAPTER VII 

CONCLUSIONS AND RECOMMENDATIONS 



7.1 CONCLUSIONS AND RECOMMENDATIONS. 


{ 


The limitation of the SIR-C/X-SAR experiment is inadequate SCR at low rain 
rates. At moderate rain rates, for off vertical look angles, there is adequate SCR if 
the angle of incidence is high (0 > 60°). At vertical due to the use of non-ideal 
weighting (raised cos 2 ) of the chirp pulse causes range sidelobes that couple 
enormous amounts of clutter to the rain signal. However, if the experiment is over 
a dense forest, due to the smaller magnitude of the forest clutter, moderate to low rain 
rates may be detected ( > 5 mm/hr). The altitude of the rain is also important; 
throughout the study it is assumed that the rain has an altitude of 3 km. 
Measurements with rain having smaller thickness are more affected at vertical. To 
overcome the clutter problem induced due to the use of the raised cos 2 pulse, we 
recommend the use of a pulse that has very low and rapidly falling range sidelobes. 

Another suggestion is to estimate the clutter and subtract the clutter from the 
combination of the signal and clutter. This is possible due to the availability of 
independent samples. The number of independent samples required to estimate the 
clutter to within 20 % is very large. However, this number can be obtained by 
degrading the resolution in both the range and the azimuth and by using multiple 
looks. Nevertheless, the resolution cannot be increased to an arbitrarily large value. 
As the resolution increases so does the size of the cell volume. Since the observable 
rain volume is fixed, the number of cells available for processing decreases. Further, 
as the cell size is increased, the probability of a cell being partially rain filled 
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increases. Partially rain filled cells have smaller returned signal power, which 
decreases the SCR. Thus, it is concluded that the minimum detectable rain rate is 
fixed to be > 10 mm at 60°. But for vertical angle of incidence rain rates > 5 mm 
can be detected. 

In obtaining an estimate of the clutter, which involves taking a separate 
measurements in absence of rain, we assume that surface echo modification due to 
impinging rain drops is negligible. Further we assume that the ocean surface wind 
speeds are the same in the rain volume as they are outside it. Sometimes these 
assumptions, as in the case of convective storms, are not true and these are the 
limitations of this approach. To overcome this limitation it is suggested that the 
estimation of the clutter be made in the same rain volume but with a fully focussed 
synthetic aperture. When full focussing is used, because of the narrow Doppler 
bandwidth used for the surface echo, the interference due to rain is minimal. To 
illustrate this consider Fig. 7.1(a), the figure shows a rain drop’s resultant velocity 
vector. Using some typical numbers to obtain the magnitude of the resultant the 
Doppler spectrum is shown in Fig. 7.1(b). Clearly the Doppler from the ground is 
wide enough apart from the Doppler due to the rain to allow separating them. 
However, the magnitude and direction of the wind determines the magnitude of the 
resultant velocity vector and hence the Doppler spectrum. Due to the Doppler shift 
caused by earth’s rotation, it may sometimes be impossible to separate the two 
Doppler components. 
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(b) The resultant Doppler. 

Fig. 7.1 The Doppler spectrum from the surface and that from the falling rain 
drops is far enough apart to separate them. 
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However, this is a better approach to estimate the clutter component than estimating 
the clutter at a location where it is not raining. 

Another approach that is feasible at vertical is try to eliminate the effect of the 
range sidelobes by transmitting a short pulse (chirped or not), but processing the 
returns without dechirping. A 3 psec pulse was used during the SIR-C mission and 
the resulting resolution is 450 meters. Though no data are available at this point, it 
is certain that range cells as high as 500 meters above the ground have little or no 
clutter associated with their signal. However, since the rain thickness is 3 to 4 km 
the number of range cells available is limited. 
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APPENDIX - A 




Table. 2.1 Laws-Parson drop-size distribution table for various precipitation rates. 


(Ulaby et, al., 1981). 
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Table. 3.1(a) 


G and H coefficients used in the SASS model. These coefficients depend 


wind direction and angle of incidence. 
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The CM0D4 Model Formulnlion and Coefficient* 


The form of ihe ECMWF CM0D4 model ii: 

= hu (l f hj co*^ ♦ hitanhhj.co*?^) 1 * 

where: 

tnd 

{ 0 if y < 0 

log v if o < y < S 
vli/3.2 if y > 5 

»„d a,0,T, k,,h end k, ik expanded u Legendre polynomial! to a total ofl8 coefficient!. 
4, U a reiiduai correction factor to k,, and it given aa a look up table aa a (unction of 
incidence angle. 

a » c ( P u + Cj F| + Cj Tj 

7 = CdP 0 + cjPi + C4P} 

0 - <1 Al + C A Pi + C * ^2 

4, = c 10 f„ + c 11 .l / + (e„F u + c„.K).r , (*) 

4j * ti(/i + cii.(l + 

4, = 0.42(1 +c„(c lf + «)(e„ + V)) 

k, = LUT(f) 

?'{z) m tank { + 2.5(* + 0.35)) — 0.6l(* + 0.35) 


where the Legendre polynomeli in * are: 


P.al P, - * r, = (3*’ - m 


with X = (» - 40)/25 


V ii the wind .peed in m.-’, * the relative wind direction in degree! and « the incidence 


angle in degree*. 


| CMOD4 Coefficient* | 


Model: 

CMOD4 

0 c y 
<1 
Cj 

2.301573 
• | .632686 
0.761210 

7 << 

Cj 
C* 

1.166619 

0.595955 

•0.293819 

0 <r 
Ca 

<» 

• 1 .0152-14 
0.342175 

• 0.500786 

by C|U 
<11 
<IJ 
<11 

0.014430 

0.002484 

0.074450 

0.004023 

h <14 

<»» 

0.148810 

0.089286 

*3 <!• 

< 1 7 
<1* 

•0.006667 

3.000000 

■10.00000 


| Reiidml Fa ctor* for CMQIX4 ' ] 


r 

6, 

r 

*r 

r 

t. 

16 

1.076 

31 

0.927 

46 

1.054 

17 

1.076 

32 

0.923 

47 

1.063 

18 

1 .075 

33 

0.930 

48 

1.052 

19 

1.072 

34 

0.937 

49 

1.047 

20 

1.069 

35 

0.944 

SO 

1.038 

21 

1.066 

36 

0.955 

51 

1.026 

22 

1.056 

37 

0.967 

52 

1.016 

23 

1.030 

36 

0.978 

53 

1.002 

24 

1.004 

39 

0.988 

54 

0.989 

25 

0.979 

40 

0.998 

55 

0.965 

20 

0.967 

41 

1.009 

56 

0.941 

27 

0.9S8 

42 

1.021 

57 

0.929 

28 

0.949 

43 

1.033 

56 

0.929 

29 

0.941 

44 

1.042 

59 

0.929 

30 

0.934 

45 

1.050 

60 

0.929 


Table. 3.1(b) Coefficients used in the CMOD-4 model. 




Table 3.2 


Wi»fi««|ih iC*»t 


Q«iAl«iy 

i*C» 

10 

l.:i 

\.u 

0 *2 



’ 20 

8.88 

8 14 

6.15 

4.44 



10 

9.02 

7.80 

5.45 

3.94 

n 


0 

8.99 

7.14 

4.75 

).45 



8 


6.48 
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Z e -R relations calculated using Mie Scattering. (Battan 1973) 
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